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^j^i Abstract. If no light Higgs boson exist, the interaction among the gauge bosons becomes strong at high energies 

r—{ ' (~ ITeV). The effects of strong electroweak symmetry breaking (SEWSB) could manifest themselves as anomalous 

couplings before they give rise to new physical states, thus measurement of all couplings and their possible deviation 
from Standard Model (SM) values could give valuable information for understanding the true nature of symmetry 

^ : 

breaking sector. Here we present a detail study of the measurement of quartic gauge couplings in weak boson scattering 
processes and a possibility for same measurement in triple weak boson production. Expected limits on the parameters 
04, 05, og, ar and oio in electroweak chiral Lagrangian are given. 



1. INTRODUCTION 

Investigation of the true mechanism that triggers electroweak symmetry breaking is one of the most intriguing 



00 

o 

• questions that are still open in the SM. If we start with minimal set of assumptions and adopt an effective Lagrangian 
approach P| then all our knowledge and ignorance about the underlying working theory is parameterized in set of 
Qh! couplings that are associated with operators making the theory finite to a given order. In the Higgs-less scenario 
O ^' anomalous quartic couplings a^, a^, ae, ay and aio of (longitudinally polarized) W and Z boson are of interest since 
1^ they give information on the underlying dynamic of the new symmetry. The set of operators that we consider 

> ■ ^4 = ^triV^VMV^^V^) , L5 = ^HV^V^)tr{V,V-) (1) 

H 

pi, - ^tr{V^V,)triTV^))tr{TVn , Lj = ^triV^V^)triTV,))triTVn , = ^(tr(TT/^))<r(rK))' (2) 

Hi 

uses the same notation as in 2], thus the results can be directly compared. Equation contains SU(2)c conserving 
operators and ones that are allowed if isospin symmetry is not conserved. Full set of possible scattering processes 
that could be used to extract the couplings is given in Tabled where we have tried to cover all of them with at least 
one decay channel (fully hadronic). Only in the multi-parameter analysis do mutual interplay of the couplings come 
in to the result, thus making the prediction more realistic. 



Table I: Sensitivity to quartic anomalous couplings for all possible scattering processes 
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Table II: Generated processes and cross sections of signal and background for ^/s — ITeV, polarization 80% left for electron 
and 40% right for positron beam 
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2. VECTOR BOSON SCATTERING 

2.1 . General Layout 

We assume a center of mass energy of ITeV and a total luminosity of lOOOfb"^ in e^e~ and 350fb~^ in e~e~ 
mode. Beam polarization of 80% for electrons and 40% for positrons is also assumed. The six fermion processes 
under study correspond to the scattering of longitudinal weak bosonos. Since triple weak boson production is also 
sensitive to quartic anomalous couplings ( ZZ or W^W~ with neutrinos of second and third generation as well 
as a part of VfAJ^WW (ZZ\ ev^W Z and e~^e~W~^W~ final states ) there is no distinct separation of signal and 
background. Signal processes in separate analysis are thus affected by all other signal processes as well as pure 
background. In comparison to the previous study Q single weak boson production was included in background 
for completeness and in order to get closer to the experimental conditions initial state radiation was taken into 
account when generating events. For the generation of tt events Pythia was used, for all other processes the full 
six fermion generator WHIZARD Q was used. No flavor summation was done since all possible quark final states 
were generated. Hadronisation was done with Pythia. The SIMDET |^ program was used to producethe detector 
response of a possible ILC detector. Table ^ contains a summary of all generated processes used for analysis and 
corresponding cross sections. For pure background processes a full lab~^ sample was generated, all signal processes 
were generated with higher statistics. Single weak boson processes and qq events were generated with an additional 
cut on M(g5)>130GeV to reduce number of generated events. 

The observables sensitive to the quartic couplings are the total cross section (either reduction or increase depending 
on the interference term in the amplitude and the point in parameter space), and modification of the differential 
corross section distributions over polar angle as well over decay angle. This is not a full set of observables but some 
sensitive event variables, for example transverse momentum, cannot be used since contribution of longitudinally 
polarized weak bosons is dropping faster then for transversally polarized wak bosons with increasing transverse 
momentum and a transverse momentum cut is an unavoidable tool to suppress background in analysis. 

2.2. Event selection 

Event selection was done using a cut based approach similar to previous analysis The general steps in the 
analysis were the use of final state e~(e+) to tag background (signal in eVeWZ case), a cut on transverse momentum, 
and missing mass and energy. Realistic ZVTOP b-tagging ^ was used when possible to enhance signal to background 
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Figure 1: The influence of ISR on result for 
W+W- scattering at SOOGeV, 68% confidence 
level contour, full line with ISR, dotted no ISR 



Figure 2: Reachable sensitivity using only the e"'"e~ mode (full 
line) and combination of e"'"e~ and e~e~ equally sharing time 
(dotted line) at ITeV, 68% confidence level contour 



separation. Finally cuts around the nominal masses of weak boson were used to accept only well reconstructed events. 
The effects of the introduction of ISR were investigated at center of mass energy of SOOGeV and W^W~ scattering 
as an example. The signal event sample was generated with and without ISR and the fit results were compared 
without taking background into consideration. The difference in the fit result, Fig.^show, as expected, deterioration 
of the result to a small extent. The significant effect of ISR is on the smearing of kinematic distributions and making 
harder signal background separation. 

2.3. Fit method and results 

Extraction of quartic gauge couplings from reconstructed kinematic variables was done with a binned likelihood 
fit. For each signal process statistics much larger than the nominal one (1000fb~^ for e+e") were generated and 
passed through the detector simulation. Each event is described by reconstructing four kinematic variables - event 
mass, absolute value of production angle cosine and absolute values of decay angle cosines of each reconstructed 
weak boson. The absolute value of the production and decay angles were used since there is no possibility to resolve 
quark antiquark and W'^W~ ambiguities. Matrix element calculation from WHIZARD was used to obtain weights 
to reweight the event as a function of quartic gauge couplings. Each Monte Carlo SM event is weighted by: 

R(ai,aj) = 1 + Aui + Bai'^ + Cuj + Da^"^ + Ea^aj (3) 

Function R{ai, aj) describes the quadratic dependence of the differential cross-section on the couplings. It is obtained 
in the following way: using the generated SM events (aj = 0, z = 4, 5, 6, 7, 10) we recalculated the matrix elements 
of the event for a set of five different points in ai,aj space and solve a set of linear equations for A,B,C,D and E. 
Due to the linear combination in which couplings contribute in amplitude 0, in any case five points are enough 
to determine the constants for weighting function. Choice of the points varied from process to process in order to 
fulfill the following conditions: distance of the point from SM value should be large enough not to come in numerical 
instabilities problem when solving equations and at the same time small enough not to come in to the region were 
phase space po|)ulation would be significantly different from the SM. Four dimensional event distributions are fitted 
with MINUIT 8\ maximizing the likelihood as a function of ai,aj taking the SM Monte Carlo sample as "data". 

L(ap,a,) = - E iV^*^(z,j, fc,01n(iV"-"Kz,j, fc,0)+ E iV"-"^ (z, j. A:, (4) 

i,j,k,l 'ij,k,l 

where i runs over the reconstructed event energy, j over the production angle, k and I over the decay 
a,ngles,N^^{i,j,k,l) are the "data" which coresspond to the SM Monte Carlo sample and N°''"°'''{i,j,k,l) is the 
sum of same SM events in the bin each weighted by R{ap, aq). Pure background events have R{ap, Oq) = 1, and for 
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Figure 3: Expected sensitivity (combined fit to all sensitive processes) to quartic anomalous couplings form lOOOfb ^ e^e 
sample full line (inner one) represents 68% and dotted (outer one) 90% confidence level a) conserved SU(2)c case b) broken 
SU(2)c case 



background coming from other sensitive processes the proper weight is taken into account. After separate analysis 
of each process from Table U a combined fit was done. Small fraction of double counted events that remains after 
single process analysis was uniquely assigned to one or another set according to the distance from the nominal mass 
of the weak boson pair ( for example WW or ZZ ). The analysis was primarily focused on e~^e^ mode since we 
expect significantly larger integrated luminosity in this mode. Effects of possible e^e" option were considered in 
following way. If available time for running the collider at a given center of mass energy is fixed and e~e~ option 
exists we can divide running time between e^e^ and e^e^ assuming that the ratio of their integrated luminosities 
is 3:1 and then do the fit to the whole data sample. A combined fit for 2ab~l e^e~ was done and compared with 
the combined result from lab^l e+e^ together with a 350fb^l e~e~ sample. The confidence level contours in Fig.|21 
show negligible difference in reachable sensitivity in these two cases. 



Table III: The expected sensitivity from lOOOfb" e"'"e~ 
sample at ITeV in SU(2)c conserving case, positive and 
negative one sigma errors given separately. 
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Table IV: The expected sensitivity from 1000fb~^e'''e~ 
sample at ITeV in broken SU(2)c case, positive and neg- 
ative one sigma errors given separately. 
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Table IIIII and Table IIVI contain results for weak boson scattering assuming integrated luminosity of lOOOfb 
e+e^ mode in SU(2)c conserving case and broken SU(2)c. 



3. TRIPLE BOSON PRODUCTION 

We consider now the reactions e+e" W~^W~Z and e'^e~ — > ZZZ. On tree level the elementary process 
producing the WW Z final state is driven by 15 Feynman diagrams. Only one of the diagrams contains the quartic 
coupling and has to be extracted from the other interfering terms. Only the part containing a longitudinal gauge boson 
is expected to give a sizable signal related to electroweak symmetry breaking. For WW Z this part is substantially 
enhanced using polarized beams. We investigate several cases: i) unpolarized, ii) 80% right polarized electrons, and 
iii) 80% right polarized electrons along with 60% left polarized positrons. For ZZZ polarization is not substantial, 
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since the standard model background is much smaUer (two diagrams). Calculations are done using the Whizard 
event generator ^j. Since the gauge bosons are short lived states they decay. Presently we consider on-shell gauge 
bosons only (narrow width approximation) and hadronize the final state using PYTHIA Q.The three-boson state 
is characterized by three four-momenta and the bosonic spins. In general three momenta lead to 12 kincmatical 
variables that are reduced by four through energy momentum conservation, by three because of the on-shell condition 
mentioned before, and by two due to rotational invariance. Hence in total three independent kincmatical variables 
are left. We choose two invariant masses of the Dalitz plot, Af^^, M'^ry/ and the angle between the beam axis 
and the direction of the Z-boson. Spin of the bosons leads to additional degrees of freedom, and we may distinguish 
longitudinal [L) from transverse (T) polarization. Presently, we do not yet analyze the bosonic spins. 

The three independent kincmatical variables lead to a three dimensional histogram. If the angle 9 is not measured 
the resulting two dimensional histogram leads a Dalitz plot. We investigate the differences on the histograms as a 
function of and a^. The observable are discretize into bins and is given by 

^^^^Slzlp^ ,5) 

where tTyfe denotes the error, and k the sums over bins of M^^, M^^, and 6. 

We use the Whizard generator to produce standard model events corresponding to a luminosity of 1000 fb~^. 
The detector efficiency is simulated using the fast simulation SIMDET j^. To reconstruct WWZ, we use the decay 
WWZ — !■ 6 jets. The dominant background is due to tt — > bbWW 6 jets. In the absence of a full simulation 
at 1000 GeV we estimate from our previous studies at 500 GeV a combined effect of efficiency and purity to be 
42%. A full simulation at 1000 GeV is presently under development. The probabilities of standard model events 
are reweighted when introducing anomalous couplings a^, a^. Since the effective Lagrangian is linear in 014, any 
observable is of second order in the parameters and can be expressed by a polynomial with five parameters. The 
parameters are determined by evaluation of N^j^^° {a4 , a^) for each event and for five pairs of fixed values (0:4, 0:5). 
By inversion A^*j;?°(a4, 0^5) is know for arbitrary values of (0:4, as), viz. 

7V*h™(a4, as) = + N.^a^ + N^^^^aj + Np^.a^ + Nf^^al + Nf^^a^a^ (6) 

for each bin i,j,k. Finally we calculate and determine Aa4(a4,a5) and Aas(a4,as) for the specific values 
— 2.30 (68.3% confidence) and — 4.61 (90% confidence). Results are shown in Fig. ^ We find that the 
sensitivity drastically increases with polarization. Sensitivity can be improved even further by utilizing meaningful 
cuts, which has not been done in the present stage of analysis and by using the information of angular distribution 
of the jets that depends on the polarization stage of the final bosons. 




Figure 4: Expected sensitivity for 04 and 05 at -^s — 1000 GeV. The lines represent 90% confidence level. Luminosity 
assumption 1000 fb~^. a) e^e~ WWZ: unpolarized case dashed line , e~ right-polarized to 80% dotted line , e"*" additionally 
left-polarized to 60% full line b) e'^e~ — > ZZZ unpolarized. 
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Figure 5: Comparison of estimated sensitivities for Q4 and 0:5 at ^ — lOOOGeV from weak boson scattering (black) ad triple 
boson production (gray). Lines represent 68% (full) and 90% (deshed) confidence level contours. 



4. Conclusion 



Expected limits on the measurement of anomalous quartic couplings from all possible weak boson scattering 
processes were presented. In Fig. |S1 we make a comparison of estimated sensitivites from weak boson scattering 
processes and ongoing triple boson production analysis. With the same integrated luminosity and 80% left e~ and 
40% right e+ polarization for scattering and 80% right e+ and 60% left e~ polarization for triple producstion we 
obtain comparable results. This shows that luminosity sharing of opposite polarization can probably lead to the 
same overall accuracy for the measurement of quartic boson couplings. The same luminosity based conclusion was 
made after comparison of e^e^ and e^e^ running modes leaving the experimental physicist several ways to achieve 
the desired precision. 
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